Scenario for Spin Glass Phase with Infinitely Many States 
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A possible phase in short-range spin glasses exhibiting infinitely many equilibrium states is pro- 
posed and characterized in real space. Experimental signatures in equilibrating systems measured 
with scanning probes are discussed. Some models with correlations in their exchange interactions 
are argued to exhibit this phase. Questions are raised about more realistic models and related issues. 
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The term "spin glass" refers to both experimental sys- 
tems and theoretical models with enough randomness 
and frustration to preclude conventional magnetic order 
0. Despite thirty years of study, the nature of pos- 
sible spin-glass phases remains controversial. Even the 
most studied model, the Edwards-Anderson (EA) Ising 
model with 7i = —J2<i j> Jij^iCj with <jj = ±1 and 
short-range couplings independently drawn from a 
symmetric distribution |2j, is still not understood. The 
droplet /scaling scenario for the spin-glass phase has only 
one pair of states and no known inconsistencies 0, U 
Parisi has hypothesized the existence of infinitely many 
equilibrium states as found in mean-field spin glasses [3| . 
But in finite-dimensional systems the nature of and rela- 
tionship between such putative states remain obscure. 

In this paper, we introduce the montage phase, a con- 
crete scenario for a spin glass phase with many incon- 
gruent states, states differing by other than symmetry in 
most of the system . By characterizing this phase and 
its experimental signatures, we also generate questions 
about more general short-range spin glasses. Since real 
systems do not have the independence and special sta- 
tistical symmetry of EA models, we consider a broader 
class of spin-glass models. For illustration, we construct 
special models that exhibit infinitely many incongruent 
states at low temperatures. 

For experiments it is essential to address a fundamen- 
tal question: How can one determine — even in princi- 
ple — the number of equilibrium states in a system too 
large to equilibrate fully? An ideal experiment would 
map repeatedly the time evolution after quenches of spin 
correlations in the same large region of size I. After a 
sufficiently long waiting time, a map of local magnetiza- 
tions in the region will resemble one of the equilibrium 
states. The number of states distinguishable in the re- 
gion that can arise from random initial conditions gives 
the number of distinct such maps, M(€). If there are in- 
finitely many incongruent equilibrium states, M{€) will 
increase with I. 

The nature of ground states and their excitations usu- 
ally determines the equilibrium behavior of low temper- 
ature phases. Infinite-system ground states are config- 
urations whose energy cannot be lowered by any finite 
change. These can be constructed using sequences of 



fixed boundary conditions on a set of nested boxes, {Bl}, 
to generate sequences of finite-region ground states. Each 
such sequence that converges yields an infinite system 
ground state. At non-zero temperature, a similar proce- 
dure yields infinite system equilibrium states, pure states. 
While most sequences will not converge, any sequence has 
at least one convergent subsequence Q. Motivated by 
Newman and Stein [tj, on whose work we build, we call 
the distribution of states that appear in the ensemble of 
sequences of random boundary conditions the equilbrium 
metastate. With short-range interactions, the number of 
possible boundary conditions on a box is exponential in 
its surface area so the number of equilibrium states dis- 
tinguishable in a box of size I grows no faster than e cl 
0- 

The physics of equilibration from random initial condi- 
tions bears some resemblance to formal limits of bound- 
ary condition sequences. After a quench, the length scale, 
L w , over which the system is in local equilibrium in- 
creases with waiting time, t w . After times larger than 
t w , regions smaller than L w resemble physical equilib- 
rium states. More precisely, in scale-i? subregions with 
I -C L w , correlations averaged over times short compared 
to t w are approximated by those of one of the infinite 
system equilibrium states 0, 0, 0|. The particular 
infinite system state that appears depends upon the sur- 
roundings of the scale-Lu, region, which itself resembles 
an equilibrated finite-sized system with effective bound- 
ary conditions imposed by these surroundings. The mat- 
uration metastate gives the distribution of equilibrium 
states that appear over the ensemble of all possible initial 
conditions and ensuing long-time dynamical histories. 

In a ferromagnet, for example, from random initial 
conditions the | and J. states appear equally frequently, 
but which occurs in a given region varies stochastically 
with t w as domain walls pass through. Typical scale- 
L w , regions will contain at most a few domain walls — 
no more than boundary conditions on such regions can 
induce in equilibrium. In principle, there are also equi- 
librium states with a fixed flat infinite domain wall, but 
these arise only from special initial conditions or special 
sequences of boundary conditions and thus appear with 
zero probability in the metastates. In any case, domain- 
wall states differ from both the | and J, states only in a 
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negligible fraction of the system and thus are regionally 
congruent |(| to these basic states. 

The equilibrium and maturation metastates of random 
Ising systems can contain more than two states. A three- 
dimensional system with all ferromagnetic couplings ex- 
cept for couplings of random sign across a single plane 
through the origin provides a simple example. The cou- 
pling between half-spaces scales as \fT? in an L x L x L 
box centered at the origin so boundary conditions can 
control each half-space separately. The resulting four 
ground states are regionally congruent because almost 
every region is the same in each up to the global sym- 
metry. The metastates contain all four states with equal 
weight since the effects of random boundary conditions 
and of random couplings are comparable. The balance 
between these determines the orientations of each half 
space in any given finite system or at any given time. 

A system with many incongruent states would, in 
contrast to one with only domain wall states, exhibit 
multiple locally-distinguishable spin maps for any large 
t w El EH 113 • Can this occur in short-range Ising spin- 
glass models? 

It is instructive to show first that some finite dimen- 
sional Ising models can have infinitely many incongruent 
states. This can occur if there are many infinite clusters 
of spins that are separately controllable. 

Consider a d-dimensional lattice consisting of a stack 
of slabs of thickness s and dimension d c = d — 1. Cou- 
plings within each slab are ferromagnetic, while those 
between them — across borders of dimension d c = d — 1 
— have random signs. The energy to change the orienta- 
tion of a size £ region within an isolated slab scales as £ ° 
with 9 C = d — 2, the cluster stiffness exponent. Inter-slab 
random coupling scales as £^ c with <f> c = d c /2. An Imry- 
Ma argument |l4j implies that each slab is ferromagnetic 
provided d — 2 > (d — l)/2. Boundary conditions can 
control the 0(£) slabs in box Bi independently. Thus 
there areTV(^) = exp(£log2/s) possible slab orientations, 
each corresponding to an infinite system state. Two ran- 
domly chosen ground states are incongruent since most 
sites are near a slab with different relative orientation 
in the two states. Isotropic versions of this model can 
be made from interpenetrating arrays of randomly cou- 
pled ci c -dimensional ferromagnetic slabs perpendicular to 
each set of d — d c directions. For d c > 2, these will have 
exponentially many incongruent equilibrium states with 
log[Af(£)] = 0(£ d ~ dc ). Although they themselves are not 
spin glasses because both their couplings and spatially 
averaged spin correlations exhibit long-range order, the 
slab models suggest an analogous spin-glass phase. 

The following fundamental property characterizes a 
putative spin-glass phase with infinitely many incongru- 
ent states which we call a montage phase: Repeated 
quenches from high temperature reveal clusters of spins 
with the same relative orientation at long times after each 
quench. The maturation metastate characterizes the set 



of such clusters of spins with similar intra-cluster ori- 
entation together with the variation of relative cluster 
orientations with wait time and from quench- to- quench. 

To exhibit a montage phase, the system must be com- 
posed of infinitely many infinite clusters of spins, each 
separately controllable by boundary conditions. Differ- 
ent states will be incongruent only if most spins are close 
to cluster boundaries, thus (i) together the clusters oc- 
cupy a finite fraction of space and (ii) their boundary 
and full dimensions are equal, d c = d c . Clusters will be 
independently controllable if (iii) in isolation each has a 
unique pair of ground states, and (iv) they interact with 
each other weakly enough. The number of distinguish- 
able states will grow with scale as Af(£) ~ exp(c£ d ~ d °). 

As with the slab models, conditions (iii) and (iv) hold 
only if the minimal energy cost to invert the spins in a 
typical scaled section of an isolated cluster is much larger 
than the typical interaction energy of this section with 
spins in all other clusters. If the stiffness scales as £ c 
and the interactions as , this requires 

0c>4>c- (1) 

The stiffness exponent, 8 C , depends on cluster topology 
and on the range and degree of frustration of the intra- 
cluster couplings. For any model with nearest-neighbor 
or sufficiently short-range interactions, d c > d— 1 so that 
clusters span the system and 

&c<d c - 1. (2) 

More precisely, 9 C < d c , the cutting dimension of a clus- 
ter: a minimal cut to sever a typical scale L-section from 
the rest of the cluster breaks of order L d ^ bonds. In gen- 
eral, dc < d c — 1 , the upper bound obtaining only if the 
clusters are very well connected. The stiffness exponent 
will be equal to d c only if bonds within a cluster are 
sufficiently unfrustrated. 

The interaction exponent <j) c depends on inter-cluster 
couplings. For independently-random couplings between 
clusters when short-range couplings dominate, 

4>c = d c /2. (3) 

Anti-correlations in the inter-cluster couplings that com- 
pensate for correlations among intra-cluster couplings 
could reduce 4> c . 

The montage phase is stable to a small uniform (or 
random) field provided that 8 C > d c /2. This follows 
by analogy to random-field ferromagnets. In a magnetic 
field the orientation of each cluster in a finite size system 
with random boundary conditions will be determined by 
the competition between the field and the inter-cluster 
energies. 

Can this spin glass phase with infinitely many states 
exist in EA models? Can it exist in other short-range 
Ising spin-glass models that are difficult to distinguish 
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from those with independently random couplings? We 
first address the latter question by constructing exam- 
ples. 

If interactions within clusters decay as a power of sep- 
aration, 

/. . - ^11 (A) 

with the Kij of order 1, cluster stiffness can be substan- 
tial even for sparse clusters and is largest for unfrustrated 
couplings within clusters. The energy to flip a typical 
scale- L region of an isolated cluster from its ground state 
orientation determines C . Since flipping a typical spin 
in an isolated unfrustrated cluster costs L da /L d+i ', sum- 
ming over pairs of spins yields 

6 C = 2d c - d - C (5) 

for C < 2d c - d. 

If the strength of inter-cluster interactions also decays 
with power law Eq. 0] but with independently random 
signs, K^ teT = ±e with equal probability, the inter- 
cluster part of the interaction per spin will be finite as 
long as C > —d/2. Then <f> c = d c /2 as in the short-range 
case. Thus such a power-law model exhibits exponen- 
tially many incongruent states if 3d c /2 > d + ( > d/2, 
which is satisfiable if d c > d/3. 

When can space-filling sets of fractal clusters exist? 
They certainly can be constructed hierarchically. For 
example, divide a d-dimensional (hyper)cube into 2 qd 
smaller cubes and label half of these A and half B. Re- 
peat this process inductively, subdividing each subcube 
and appending additional labels A and B. Each string 
of As and Bs will then label a fractal cluster of dimen- 
sion d c = d — 1/q. Furthermore, provided that d > 3 
and q is sufficiently large, this can be done with some 
randomness in the labeling but so that nevertheless each 
cluster is connected and system spanning with boundary 
dimension d c = d c and cutting dimension dc > d c /2. 

These hierarchically constructed clusters with nearest 
neighbor interactions that are unfrustrated within clus- 
ters and random in sign between them are separately 
controllable. In this case, 9 C — d) c > d c /2 = cj> c so that 
nearest neighbor couplings are sufficient to generate ex- 
ponentially many states in three (or more) dimensions. 
A set of signs labels each state- one sign for every infinite 
string of As and Bs . 

More realistic models require space-filling sets of ap- 
propriate fractal clusters that are statistically transla- 
tionally invariant. Do these exist? Following Newman 
and Stein fl(| , one possibility is built from invasion per- 
colation clusters |2£|. In dimensions d > 8, the number of 
disjoint invasion clusters is believed to grow as a power of 
the scale. Similarly, loop-erased random walks in d > 4 
form sets of clusters with d c > 2 that are collectively 



space filling [Tjj. With power law interactions with spe- 
cial correlations, such clusters can be separately control- 
lable. With short-range interactions, however, it is not 
known whether they will be, even in high dimensions. 

The special interactions in the above power-law ex- 
ample are physically unrealistic. It seems unlikely that 
conventional EA models exhibit separately controllable 
clusters because of the apparent need for many weakly 
frustrated sets. In real spin glasses, however, couplings 
are neither symmetrically distributed nor independent 
of other nearby couplings, and can be long range: e.g. 
RKKY (~ 1/r 3 ). Can clusters satisfying conditions (i)- 
(iv) arise spontaneously in systems with interactions that 
have translationally invariant correlations among them? 
If so, how short-ranged can the interactions and correla- 
tions be? In particular, can the correlations among the 
couplings come from an effective Boltzmann-likc proba- 
bility measure with only short range interactions in the 
effective Hamiltonian? 

XY models with infinite uniaxial anisotropy suggest 
some reason for optimism: they both illustrate a pos- 
sible mechanism for generating many states as well as 
provide an example in which controllable fractal clusters 
arise spontaneously. In dimensions greater than four, XY 
models with random uniaxial anisotropy are believed to 
exhibit ferromagnetic order [lsj . When anisotropy is in- 
finite, the random axis at site i, a», dictates the direction 
of the XY spin: Si = (Ji&i with cr^ = ±1. The result- 
ing Hamiltonian is equivalent to an Ising Hamiltonian 
with Jij In high dimensions, this model should 

exhibit ferromagnetic order in the XY variables provided 
that the coordination number is sufficiently high or corre- 
lation in random axis directions is sufficient [19( . Conse- 
quently the direction of the spontaneous magnetization, 
fx, should parametrize an infinite family of states — with 
<7jaj • fj, > at most sites — in a way which would be 
hidden in the Ising variables. 

Although this may be a qualitatively distinct scenario 
for many states, the number of which probably would 
grow only as a power of length scale some separately con- 
trollable fractal clusters do exist in such models. All spins 
in a given finite box flip when fixed boundary conditions 
inducing a spontaneous magnetization /i are rotated by 
7T. There must be some much smaller change in bound- 
ary conditions inducing a new state, // that is as close as 
possible to /i while still differing in the deep interior of 
the box. Reversal of a connected fractal set of spins that 
boundary conditions can control relates states fi and //. 
But interactions between such sets are too strong for this 
system to be in the montage phase. We expect that the 
number of controllable fractal sets only grows as a power 
of system size and that at most only a few of them can 
be controlled simultaneously. 

We now turn to experimentally observable conse- 
quences of the montage phase scenario. Maps of local 
magnetizations (or other local properties) in the same 
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region after repeated quenches — to the same conditions 
followed by waits for time t w give the spectrum of 
states in the maturation metastate as well as their local 
structure. Soon after a quench from high temperature, 
only nearby spins are correlated and there will be e l l v 
distinct possible sets of spin orientations. With increas- 
ing t w , V will increase, and the number of observed maps 
will decrease. If there are only two states, after a suffi- 
ciently long time such maps will coincide in the observed 
region up to a global spin inversion. But in the montage 
phase, even after extremely long times, the number of 
distinct states observed in a scaled region will grow with 
t as 

M{1) ~ exp(r) (6) 

presumably with a = d — d c . Though estimating the 
number of realizable states may be impracticable, exper- 
iments should be able to determine whether a > 0. 

Correlations between runs will carry information about 
the structure of physical states. Long-distance correla- 
tions between pairs of spins on the same cluster, (SiSj), 
averaged over a long time (but short compared to t w ), 
will be reproducible from run-to-run whereas the sign 
of correlations between spins on different clusters will 
vary. So intra-cluster correlations will dominate the 
mean-square correlations obtained by first averaging over 
many runs to obtain Y xy = (S x S y ) and then averaging 
T^y over many pairs of spins distance i apart. These cor- 
relations decay as U p to the length scale L w (t w ) 
on which equilibrium is established. Decay is exponential 
for I 3> L w . Spins in the same cluster have T xy = 0(1). 

There are subtleties in determining whether local equi- 
librium has been reached on a sufficiently large scale L w 
relative to I. First, pairs of spins on neighboring clusters 
may appear to be on the same cluster if these clusters 
tend to have the same mutual orientations at early times. 
Second, two sites on the same cluster may appear as if 
they are on distinct clusters if they are linked primarily 
via other spins that that extend well beyond L w . Both 
of these could affect the inferred a as well as the decay 
of averaged correlations. 

We have proposed a scenario for a spin-glass phase 
with exponentially many states that would have defi- 
nite experimental signatures. Such a phase can occur 
in spin-glass-like models with correlations between cou- 
plings that are special but would be hard to detect in 
a sample. Whether it can occur in less artificial models 
is unclear. But the scenario, if it does occur, gives rise 
to many questions. What are appropriate measures for 
the necessary degree of frustration in couplings? How 
do states change with temperature and magnetic field? 
What are the natures of phase transitions into such a 
phase? 
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